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INTRODUCTION 


This  report  presents  measurements  of  temperature  in  the  upper  ocean 
obtained  by  use  of  a  towed  thermistor  chain.  The  measurements  were  taken 
as  part  of  the  Joint  Air-Sea  Interaction  (JASIN)  Experiment  conducted 
during  the  summer  of  1978.  A  summary  of  the  scientific  and  operational 
plans  for  JASIN  has  been  given  by  Pollard  (1978).  More  detailed  accounts 
are  given  in  documents  published  by  the  Royal  Society  (1977,  1978). 


INSTRUMENTATION 

The  thermistor  chain  consisted  of  sensors,  electrical  conductors, 
plastic  fairing,  a  strain  member  and  a  450  kg  lead-filled  depressor.  The 
thermistors  were  manufactured  by  Thermometries  (Model  P-85).  They  were 
molded  into  sections  of  fairing  together  with  bridge/amplifiers.  Power 
and  signals  were  transmitted  by  electrical  conductors  running  through  the 
tail  sections  of  the  fairing.  The  thermistors  were  spaced  at  intervals  of 
2  m  over  a  section  of  chain  50  m  in  length.  Four  pressure  sensors  were 
also  installed  on  the  chain  at  15  m  intervals.  The  pressure  sensors  were 
manufactured  by  Kulite  and  were  installed  in  tail  sections  together  with 
bridge/amplifiers  in  a  fashion  similar  to  the  thermistor  electronics. 

Signals  from  the  sensors  were  recorded,  processed  and  displayed  by  use  of 
a  minicomputer  system  manufactured  by  Digital  Equipment  Corporation  (PDP 
11/05).  A  more  complete  description  of  the  thermistor  chain  system  is  given 
by  Spoering  and  Paulson  (1980). 
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OBSERVATIONS 

The  thermistor  chain  was  towed  by  the  R/V  ATLANTIS  II  in  an  area  about 
400  km  northwest  of  Scotland.  The  locations  of  the  tows  are  tabulated  and 
plotted  in  the  section  entitled  Tow  Tracks.  A  summary  of  the  tow  parameters 
is  given  in  Table  1.  The  first  digit  of  the  run  number  designates  the 
number  of  the  deployment  of  the  chain.  With  the  exception  of  the  first 
deployment,  the  chain  was  always  towed  counterclockwise  around  a  square. 

The  letter  following  the  deployment  number  designates  the  side  of  the  square, 
i.e.  N  indicates  a  tow  leg  toward  the  west  on  the  north  side  of  the  box. 

The  digits  following  the  dash  in  the  run  number  designates  the  leg,  one  leg 
per  side  of  the  square  except  for  the  first  deployment.  The  first  deployment 
was  separated  into  three  legs.  The  middle  leg  contained  a  front  which  sep¬ 
arated  regions  of  nearly  constant  thermal  properties.  The  square  around  which 
tows  were  made  was  usually  15  km  on  a  side  surrounding  the  Fixed  Intensive 
Array.  However,  on  two  occasions  tows  were  conducted  in  cooperation  with 
other  ships.  The  first  occasion  was  around  a  five-km  square  containing 
the  drifting  buoy  PI  and  the  second  occasion  was  around  a  five-km  square 
containing  the  mooring  H2.  The  instrumented  section  of  the  chain  extended 
from  about  20  to  70  m  depth  in  all  tows  except  the  last  which  was  10  m 
shallower.  The  tow  speed  was  usually  about  3  m/s.  The  tow  speeds  tabulated 
in  Table  1  were  determined  from  L0RAN  C  and  satellite  navigational  fixes. 


TABLE  1.  Thermistor  chain  tows  during  JASIN-78. 
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ANALYSIS 

The  temperature  observations  were  low-pass  filtered  by  computing 
sequential  30  s  averages.  Filtering  removes  variations  caused  by  surface 
gravity  waves  and  ship  heave.  The  filtered  observations  are  shown  in  the 
section  entitled  Temperature  Cross-Sections. 

Isotherm  depths  were  determined  by  linear  interpolation  between  the 
filtered  temperature  observations.  The  depths  of  isotherms  in  intervals 
of  0.2°C  are  shown  in  the  section  entitled  Isotherm  Cross-Sections. 

Spectra  of  the  shallowest  and  deepest  isotherm  on  each  leg  were 
computed  and  are  presented  in  the  section  entitled  Spectra  of  Highest  and 
Lowest  Isotherms. 

Spectra  of  the  depth  of  the  thermistor  chain  at  three  locations  on 
the  chain  are  shown  in  Figure  1.  The  spectra  were  computed  from  low-pass 
filtered  (30  s  averages)  measurements  of  pressure  during  Run  1-1.  The 
magnitude  of  vertical  oscillations  of  the  chain  is  usually  greatest  near 
the  bottom.  The  magnitude  of  the  spectra  can  be  compared  with  spectra  of 
isotherm  depth  shown  in  the  section  entitled  Isotherm  Cross-Sections.  The 
magnitude  of  the  isotherm  spectra  are  usually  more  than  two  orders  of 
magnitude  greater  than  spectra  of  depth.  We  therefore  conclude  that  variations 
in  depth  of  the  chain  have  a  negligible  effect  on  spectra  of  isotherm  depth 
determined  from  the  low-pass  filtered  temperature  measurements. 

More  details  on  analysis  procedures  can  be  found  in  Spoering  and 


Paulson  (1980). 


RUN  1-1  24--RUG-78 
PRESSURE  SPECTRA 


-5  -U  -3  -2 


CPM 


Figure  1.  Spectra  of  depth  measured  at  three 
locations  on  the  towed  thermistor 
chain  during  Run  1-1.  The  mean 
depths  of  the  pressure  sensors  are 
indicated.  The  pressure  records 
were  low-pass  (30  s  averages) 
filtered  prior  to  computing  spectra 
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TOW  TRACKS 

On  the  following  pages  there  is  a  tabulation  of  positions  during  each 
of  the  tows  followed  by  plots,  one  for  each  run.  The  tabulated  positions 
are  plotted  (x)  on  the  plots.  The  positions  were  determined  by  use  of 
LORAN  C  and  navigational  satellite.  The  local  coordinates,  x  and  y,  are 
in  kilometers  north  and  east,  respectively,  of  a  point  at  the  center  of 
the  Fixed  Intensive  Array  (59°N,  12.5°W).  The  locations  of  moorings  B1 , 
B2,  B3  and  H2  are  shown  in  several  plots. 
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RUN  1  24-AUG-78 
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X 

V 

1609 

59 
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12 

26.  09 

3.  75 
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59 

3  54 

12 
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6.  57 
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59 

3.  37 

12 
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59 
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12 
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59 
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12 
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1.  98 

13 

21.  52 

-49.  35 

3.  63 

2101 

59 

1.  88 

13 
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2  25-AUG-78 
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X 
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1203 
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12 
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40.  25 
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4.  60 

1304 
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53.  32 
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-9.  14 

-3.  12 
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55.  75 

12 
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-8.  43 
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1447 
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12 

24.  68 
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58 
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12 
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12 
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59 
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12 
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7.  91 
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59 

3.  46 

12 
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1741 

59 
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12 
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5.  37 
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59 
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12 
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2.  53 
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12 
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58 
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12 
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RUN  5  31--RUG-78 
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COORDINATE 
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RUN  ?  4-SEP-73 


LfiTITUDE 

LONGITUDE 

LOCRL 

COORDINRTES 

T  I  ME 

DEG 

MIN 

DEG 

MIN 

X 
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7.  41 

1200 

59 

3.  55 

12 

38.  OO 

-7.  66 
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-7.  96 

-7.  33 
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3.  93 
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7.  30 

1718 

59 

3  99 

12 

38.  45 

-8.  09 

7.  41 

1847 

58 

56  10 

12 

38.  40 

-8.  05 

-7.  24 
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TEMPERATURE  CROSS-SECTIONS 

On  the  following  pages  there  are  plots,  one  for  each  leg,  of  low-pass 
filtered  temperature  as  a  function  of  time  and  distance  along  the  leg.  The 
filtering  was  accomplished  by  computing  sequential  30  s  averages.  The  mean 
depth  of  each  sensor  is  given. 
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ISOTHERM  CROSS-SECTIONS 


On  the  following  pages  there  are  plots  of  isotherms  at  0.2°C  intervals 
which  were  obtained  by  linear  interpolation  between  the  low-pass  filtered 
temperature  observations  shown  in  the  previous  section.  Isotherms  which 
were  not  at  least  80  complete  were  not  plotted.  Isotherms  which  were 
incomplete,  but  had  no  more  than  20-  of  the  record  missing,  were  completed 
by  linear  extrapolation  from  adjacent  isotherms. 
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RUNS  (-N- 1  -  PU-2  29-HUG-78 
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SPECTRA  OF  HIGHEST  AND  LOWES.  ISOTHERMS 
Spectra  of  the  depth  of  the  shallowest  and  deepest  isotherms  from  each 
leg  are  plotted  on  the  following  pages.  The  symbol  x  designates  the 
va  1  lowest  i  .->0  then:  (hiohr-.t  tempera  lure)  and  H  designates  the  deepest 
isotherm.  The  rime- series  o-  low-pass  filtered  isotherm  depth  were 
prewhitened  by  taking  first  differences  prior  to  computing  Fourier  coef¬ 
ficients.  The  spectra  were  then  recolored  and  smoothed  by  averaging  over 
non-overlapping  frequency  bands,  five  per  decade,  equally  spaced  on  a 
logarithmic  scale. 
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The  purpose  of  this  Appendix  is  to  derive  an  expression  for  the  shape 
of  the  chain  under  tow.  The  derived  shape  serves  as  the  basis  for  inter¬ 
polating  the  depths  of  the  sensors  between  pressure  measurements  and  for 
smoothing  variations  in  the  pressure  records  caused  by  errors  in  the 
measurements . 

We  simplify  the  problem  by  considering  the  equilibrium  case  in  which 
accelerations  are  neglected  and  the  shape  is  independent  of  time.  The 
problem  then  becomes  similar  to  the  problem  of  the  "hanging  chain",  treat¬ 
ments  of  which  are  given  in  many  textbooks  (e.g.,  Sokolnikoff  and  Redheffer, 
1958,  pp.  40-42). 

A  schematic  diagram  of  the  chain  under  tow  is  shown  in  Figure  A1 .  The 
x-z  coordinate  system  has  its  origin  at  the  bottom  of  the  chain  in  the 
center  of  the  dead-weight  depressor  (bomb).  The  forces  acting  on  the  chain 
at  x  =  z  *  0  are  Gg,  the  gravitational  force  on  the  bomb;  Dg,  the  drag  force 
on  the  bomb;  T  (s  =  0),  the  upward  tension  in  the  chain,  where  s  is  the 
distance  along  the  chain  from  x  =  z  =  0.  The  forces  acting  on  an  element 
of  the  chain,  6s,  at  an  arbitrary  point  s  on  the  chain  are:  W(s)6s,  the 
gravitational  force,  where  W(s)  has  dimensions  of  force  per  unit  length; 
D(s)6s,  the  drag  force,  where  D(s)  also  has  dimensions  of  force  per  unit 
length;  T(s  +  6s/2)  the  upward  tensile  force;  and  T(s  -  <Ss/2),  the  downward 
tensile  force.  The  buoyancy  forces  are  incorporated  into  the  gravitational 
forces,  Gd  and  W.  The  tensile  force  may  be  decomposed  into  vertical  and 

D 

horizontal  components  as  shown  where  e ( s )  is  the  angle  between  the  chain 
and  the  horizontal . 

If  the  chain  is  in  equilibrium,  we  may,  as  shown  in  Figure  A1 ,  separately 
equate  the  horizontal  and  vertical  forces  acting  on  the  chain  at  any  point 
along  the  chain: 


T'~ 
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DR  +  /  D(s)  ds  =  T(s)  cos  e{s) 


GD  +  f  W(s)  ds  =  T(s)  sin  e(s) 

B  0 

We  now  assume  that  D(s)  and  W(s)  are  independent  of  s.  This  is  essentially 
true  for  W(s)  and  is  approximately  true  for  D(s)  providing  e(s)  does  not 
differ  too  far  from  90°.  The  equations  (Al)  are  then  integrated  to  obtain: 


Dg  +  sD  =  T(s)  cos  e(s)  (A‘ 

Gg  +  sW  =  T(s)  sin  e(s)  (A; 

T(s)  can  be  eliminated  from  (A2)  to  obtain: 

Gb  +  SW 
tan  6  ( s )  =  dTTTd 

D 

It  is  now  possible  to  solve  for  e(s),  x(s)  and  z(s).  e(s)  is  obtained 
directly  from  (A3)  and  x(s)  and  z(s)  are  given  by: 

S 

x(s)  =  /  cos  e( s 1 )  ds' 

0 


z(s)  =  /  sin  e (s ’ )ds ' 

0 

Consider  first  the  expression  for  z(s).  Substituting  from  (A3)  into  (A4) 


yields: 


z(s)  =  / 


(GB  +  s 'W)ds ' 


0  [(Gr  +  s*W)2  +  (Dr  +  s'D)2] 


which  can  be  rewritten 


s  (Gn  +  s'W)ds' 
z(s)  =  /  - 5 - 

0  (a  +  bs '  +  cs'2) 


f  , 
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2  2 

where  a  -  Gg  +  D 

b  ?  26gW  +  2  OgD 

2  2 

c  5  VT  +  D 


Integrating  (A5)  yields: 


Gq  bW 

ZU)  *  ^  -  2c 3/2 


[sinh'1 


✓d 


-  si nh"1  (^)]  (A7) 

l, 

+  *  [(A  +  bs  +  cs2)  -  /a] 

where 

d  =  4ac  -  b2 

An  expression  for  x(s)  can  be  obtained  by  a  procedure  similar  to  that 
used  to  obtain  (A7).  From  (A3)  and  (A4)  we  obtain 

s  (0B  +  s'P)ds' 

x{s)  =  /  T4  (A8) 

0  (a  =  b  s '  +  cs '  ) 

where  a,  b  and  c  have  the  definitions  given  in  (A6).  If  one  replaces  Gg 
and  W  in  (A5)  with  Dg  and  D  respectively,  the  equation  becomes  identical 
to  (A8).  The  integration  of  (A8)  therefore  yields  an  expression  of  the 
same  form  as  (A7)  with  Gg  and  W  replaced  by  Dg  and  D: 

x ( s )  =  l-p  -  — j-rp )  [sinh-1  /_2cs_+_b\ 
vc  2c  7  yj 

-  sinh"1  (-^-)] 

/d 

+  ^  [(a  +  bs  +  cs2)  -  </a  ] 


(A9) 
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The  drag  forces,  Dg  and  D,  are  assumed  to  be  proportional  to  pi) 
where  p  is  the  density  of  the  water  and  U  is  the  tow  velocity: 


DB  '  CB  AB  pU 


(A10) 


D  =  C  A  PU2 

where  Cg  and  C  are  the  drag  coefficients  for  the  depressor  and  the  faired 
chain,  respectively.  Ag  and  A  are  the  corresponding  cross-sectional  areas 
where  A  has  dimensions  of  area  per  unit  length. 

The  drag  coefficient  for  the  depressor  was  assumed  equal  to  one.  This 
value  can  be  compared  with  the  value  for  a  circular  disk  equal  0.55  in  the 
range  of  Reynolds  numbers  likely  to  be  encountered  (Batchelor,  1967,  p.  341). 

A  larger  value  was  assumed  for  the  depressor  because,  although  it  is  stream¬ 
lined,  it  has  protruding  fins,  cable  attachments  and  handling  rings  and  may 
at  times  present  a  larger  cross-sectional  area  to  the  flow  than  assumed.  It 
turns  out  that  the  shape  of  the  chain  is  not  critically  dependent  on  the 
assumed  value  of  Cg.  If  one  decreases  Cg  by  50%  at  a  tow  speed  of  3  m/s,  the 
computed  change  in  depth  95  m  below  the  surface  is  only  0.3  m.  The  Reynolds 

5 

number  for  the  bomb  at  a  tow  speed  of  3  m/s  is  6.7x10  . 

The  drag  coefficient  for  the  faired  chain  was  taken  equal  to  0.13  as 
suggested  by  the  manufacturer.  For  comparison,  bluff  bodies  have  drag 
coefficients  about  equal  0.5.  The  value  C  =  0.13  yielded  chain  shapes  in 
good  agreement  with  the  pressure  measurements  on  the  chain.  The  uncertainties 
of  the  calculated  depths  of  the  thermistors  are  estimated  not  to  exceed 

4 

±  1  m.  The  Reynolds  number  for  the  fairing  is  6x10  at  a  tow  speed  of  3  m/s. 

Chain  shapes  calculated  by  use  of  (A7)  and  (A9)  are  shown  in  Figure 
A2.  For  a  given  tow  speed,  the  shape  depends  only  the  the  distance  along 
the  chain  from  the  depressor.  The  parameters  used  to  obtain  the  curves  shown 
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in  Figure  A2  were: 

Gg  =  4,120  Newtons 
W  =7  Newtons/m 

2  _2  (All) 

CgAgp  =  12.3  Newtons/m  s 

3  -2 

CAp  =  2.6  Newtons/m  s 

The  above  parameters  are  the  best  estimates  for  the  chain. 

An  additional  model  for  the  shape  of  the  chain  was  developed  to  compare 
two  different  assumptions  for  the  drag  force  on  the  chain.  In  the  previous 
development  we  assumed  the  drag,  D,  proportional  to  s,  the  distance  along 
the  chain  from  the  depressor.  Alternatively,  we  assume  that  D  is  proportional 
to  z,  the  vertical  space  coordinate.  The  equation  corresponding  to  (A3)  is: 


dz  Gr  +  sW 

37  =  tane(s>  '  d-T7d 

D 


(A12) 


To  ease  the  integration  of  (A12)  we  assume  that  W  is  negligible,  from  which 
follows: 

x  =  (DB  +  !  z)  it  (A13) 

8 

We  may  also  derive  s(z)  as  follows: 


dz 

ds 


sin  9 ( s ) 


where 


1 

sin  e(s) 


dz' 


s  = 


iV  *  idb 


dzT] 


B 


dz’ 


z  2  ‘i 

s  =  f  (e  +  fz'  +  gz'  )  dz' 

0 


> 
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e  l  +  — y  ,  f  - 2  »  9 

h  % 


Integrating  we  obtain: 


(2gz  +  f)  /Z  -  f/e 

4g 


t?_l_  [sinh-'  (Sup £)  -  sinh-'(-|)] 

where 

Z  =  e  +  fz  +  gz^ 
q  =  4  eg  -  f^ 
k  « 

q 

The  two  models  are  compared  in  Figures  A3  and  A4.  In  Figure  A3  the  shape 
of  the  chain  given  by  the  second  model  (D  “  z)  is  plotted  for  various  tow 
speeds.  In  Figure  A4  the  shape  of  the  chain  given  by  the  first  model  (D  *  s) 
is  plotted  for  the  same  tow  speeds.  The  parameters  used  in  constructing  the 
cruves  shown  in  Figures  A3  and  A4  are  identical  to  those  used  in  generating 
the  curves  shown  in  Figure  A2  except  that  W,  the  gravitational  minus  buoyancy 
force,  is  assumed  zero  in  Figures  A3  and  A4.  The  difference  between  the  two 
models  for  a  tow  speed  of  3  m/s  amounts  to  only  0.3  m  difference  in  z  at 
s  =  70  m.  We  therefore  conclude  that  the  calculated  shape  of  the  chain  is 
not  critically  dependent  on  the  model  for  the  tow  speeds  employed  in  JASIN. 

For  faster  +ow  speeds  the  differences  may  become  significant.  For  example, 
at  a  tow  speed  of  5  m/s  the  difference  in  z  when  s  =  100  m  is  4  m. 


Figure  A3.  The  shape  of  a  towed  thermistor  chain,  l 00  m  in  length, 
for  tow  speeds  ranging  from  1  to  6  m/s.  The  chain  is 
assumed  to  be  weightless  in  water  and  the  drag  force  is 
assumed  proportional  to  z. 


